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ABSTRACT 

Aims. In this paper we re-visit the observational relation between X-ray luminosity and temperature for high-z galaxy clusters and 
compare it with the local Lx-T and with theoretical models. 

Methods. To these ends we use a sample of 17 clusters extracted from the Chandra archive supplemented with additional clusters 
from the literature, either observed by Chandra or XMM-Newton, to form a final sample of 39 high redshift (0.25 < z <1.3) objects. 
Different statistical approaches are adopted to analyze the Lx-T relation. 

Results. The slope of the Lx-T relation of high redshift clusters is steeper than expected from the self-similar model predictions 
and steeper, even though still compatible within the errors, than the local Lx-T slope. The distant cluster Lx-T relation shows a 
significant evolution with respect to the local Universe: high-z clusters are more luminous than the local ones by a factor w 2 at any 
given temperature. The evolution with redshift of the L^-T relation cannot be described by a single power law nor by the evolution 
predicted by the self-similar model. 

Conclusions. We find a strong evolution, similar or stronger than the self-similar model, from z = to z < 0.3 followed by a much 
weaker, if any, evolution at higher redshifts. The weaker evolution is compatible with non-gravitational models of structure formation. 
According to us a statistically significant sample of nearby clusters (z<0.25) should be observed with the current available X-ray 
telescopes to completely exclude observational effects due to different generation detectors and to understand this novel result. 

Key words, galaxies: clusters: general - galaxies: high redshift - cosmology: observation - intergalactic medium - X-rays: galaxies: 
clusters 
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O ' 1. Introduction 

Clusters of galaxies represent the largest scale of fully collapsed structures in the Universe and thus offer a unique insight into the 
sjj^ formation of structures and into the parameters governing their evolution. As powerful X-ray emitters (Lx - 10**^ - lO'^^erg s"') 
galaxy clusters can be easily traced up to high redshift, hence they are a very important observational tool for cosmologists. One 
problem that must be examined is how to relate the observable quantities (e.g., cluster X-ray gas temperature or luminosity) to 
|h the quantity predicted by the cosmological studies (usually mass). A detailed and accurate knowledge of the relations between 
observable bulk features is necessary. For instance, the relation between the X-ray luminosity and temperature of the intra-cluster 
medium (ICM), that is analyzed in detail in this paper, provides a nieans to convert an easily observed X-ray luminosity function 
into a more cosmolo gically useful mass function (e.g.. lEvrard et alJll996l: iReiprich & Bohringerll2002l; [Evrard et al.l 120021 : lAllenI 
l2003tlBorganill2006h . 

Under the assumption that the baryons follow the total matter distribution (dominated by dark matter) iKaiseii ( Il986l) constructed the 
simple self-similar model of cluster formation where the ICM is driven only by gravitational processes like shocks and adiabatic 
compression. As long as the baryons are distributed in the same way as the total mass, each X-ray observable scales like some 
power of t he mass. The self-simil ar model, supported by N-body hydrodynamical simulations which consider only gravitational 
processes jBrvan & Normarjl 19981) . predicts the relation between luminosity and temperature to be Lx T^. 

The observational tight corre lation between luminosity and temperature found for low redshift cluster samples (e.g iMarkevitchl 
ll998tlArnaud & Evrardlll999h indicates a similar formation history for all clusters and is qualitatively similar to the predictions of 
the self-similar models of the cluster formation. However the observed slope of the L^-T relation is found to be much steeper (^ 3) 
than expected from the self-similar model. Such deviation from self-similarity is taken as evidence that the simple gravitational 
collapse is not the only process that governs the heating of baryons, and that non-gravitational processes occurring before or during 
the cluster formation contribute to the clusters energy budget in a non negligible way. To explain the departure from self-similarity 
an increase in the central entropy is required (Evrard & Henrv 1991; Kaiser 1991; Ponman et al. 1999; Tozzi & Norman 200ll 
e.g.). This extra entropy makes the gas harder to compress in the cluster core, particularly in the shallower potential wells of low 
temperature clusters, whose thermal bremsstrahlung emission is reduced as compared to the self-similar model prediction. The 
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physical processes that can rise the entropy of the gas involve hea t ing by supernovae a nd by active galactic nuclei (AGNs)(e.g. 
iBalogh et al.l ll999'; 'Bower et al."200r; Brighenti & Mathews '2001'; 'Borga ni et alj |2002|) . or the removal of low-entro py gas via 
coohng (e.g.jB ryan 2000; Voit & Bryan 2001; Wu & Xue 2002) or a combination of heating and radiative cooling (e.g. lVoit et alJ 
I2OO2I; llbrnatore et al.ll2003i: iBorgani et al.ll2004tlMuanwong et alj|2006l) . Therefore an accurate analysis of the Lx-T relationship 
is also important to investigate the physics behind galaxy cluster formation. 

Further information on the formation and evolution of structure in the Universe are given by the study of the statistical properties of 
high redshift cluster samples with respect to the low redshift samples. Any observed change in the Lx-T scaling law at high redshift 
with respect to low redshift allows in principle to verify the predicti ons of the cluster self-similar model or to constrain different 
scenarios of heating schemes and cooling efficiency (e.g. lVoitll2005al) . 

With the main goal to revisit the L^-T relationship for high-z clusters we used seventeen clusters in the range 0.25 < z < 1.01 
selected from the Chandra archive. We analyzed the data as described in details in a companion paper jBranchesi et al.ll2007al 
Paper I). In order to increase the statistics the sample was supplemented with additional clusters from the hterature with redshift up 
to z~1.3, either observed by Chandra or XMM-Newton. The final sample consists of 39 objects. The observational results from the 
Lx-T analysis of the 39 distant clusters are compared with those obtained from the analysis of the L^-T relation of low-z clusters 
(lMarkevitchll998l;lArnaud & Evrardll 19991) and with theoretical models for structure formation. 

All the uncertainties in the paper are at the 1 cr confidence level, unless otherwise noted. We use a ACDM cosmology with Hq = 
70 km s ' Mpc ' and = 1 - = 0.3. 



2. Cluster Scaling Relations 

2. 1 . The self-similar model 

The self-similar model detailed by iKaiserl (Il986h is the most natural expectation in a picture where clusters form via the collapse 
of the most dense regions and cluster baryons are heated only by gravitational processes (adiabatic compression and shock heating) 
during the collapse. In this scenario it is assumed that a constant fraction of the baryonic content of the Universe ends up in the 
intracluster gas in hydrostatic equilibrium within the dominant dark matter, and emits by thermal bremsstrahlung. A state of virial 
equilibrium is achieved within a virial radius Ry(z) which includes a mean density which is a factor Av(z) above the critical density 
of the Universe, p^,; - 3H^ /(EttG), at that redshift. A cosmology factor E. is used hereafter to indicate the evolution of the Hubble 
constant at redshift z with respect to z=0 for a flat cosmology with matter density Q.^: 

E, = HJH^, = +z? + l- Qm]'^' (1) 

In a flat A = cosmology E^-{\+ z)'-^. 
The density contrast is given by 

Av(z) =187^2 + 82(Q, - 1) - 39(Q, - 1)^ (Bryan & Norman 1998) (2) 
with Q; = Qm(l + zy' IEI. In a flat A = cosmology Ay = 187r^ and it is independent of redshift. 

In this simple scenario, called the self-similar scenario, the properties of clusters of different masses and at different redshifts are 
related to one-another according to simple scaling laws. The expected properties of clusters at high redshift are identical to those of 
low-redshift clusters, apart from scaling factors reflecting the increase of the mean density of the Universe with redshift. The self- 
similar scenario predi cts how the ICM physical properties, like luminosity and temperature, are related (e.g. Kaiser 1986, 199l|; 
lEvrard&Hen"rvlll991l) . 

Under the assumption that the ICM is spherically symmetric, for a y6-model dCavaliere & Fusco-Femianoll976h the total mass within 
a radius R is given by the following equation (see e.g. .Voit.2005a) : 

3/3kBT(°K)Rc (R/Rc)^ 

M,ot(< R) = ^ ' , ' ' ,2 (3) 
Gyump 1 + (R/Rc)^ 

where the core radius Rc and p parametrize the gas density profile, yu is the mean molecular weight (yu=0.6 for a primordial com- 
position with a 76% fraction contributed by hydrogen) and mp the proton mass. The mass considered in the self-similar model, 
expressed as a function of the radius 7?a. which encloses a mean density which is a factor A, above the critical density of the 
Universe, is given by: 

M,ot(< RaJ = ^PcAz^l (4) 
and: 

A, = A(0)Av(z)/Av(0) (5) 
where A(0) is the density contrast at z=0. From Eq.[3]and|4]the radius corresponding to A. is calculated as follows: 

Y jj-mp HfA, 
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where /fj is defined in Eq.[T] For A, - Ay{z) one obtains the virial radius Rv(z)- These radii are oc T"^ and decrease with redshift in 
a way dependent on the assumed cosmological model. 

At very high energies the ICM behaves as a fully ionized plasma and thus it emits by thermal bremsstrahlung in the X-ray band. 
The emissivity for this thermal process at frequency v scales as: 

Ignoring the temperature dependence of the Gaunt factor g(v, T) (see, e.g.. ISpitzer II 19781) . and assuming a constant ratio of gas 
density to total mass density, the scaling behavior of a cluster X-ray bolometric luminosity can be written as: 

Uoi ^ Mpg,,T'/2 

Then one can obtain the following scaling relations between the observed bolometric luminosity Lboi and temperature of the gas T 
and the total gravitating mass: 

Uoi <x T^E.A^^ and U^i M'^^^eI^'a'J'' (9) 

The relation between the two X-ray observables, temperature and luminosity, is a useful statistical tool that enables the study of the 
physics of clusters and allows the verification of the self-similarity model expectations. 

2.2. Evolution with redshift of the Lboi - T relation 

When a high-redshift cluster sample is analyzed, it is important to define what is meant by evolution, since different definitions are 
found in the literature. The self-similar models predict changes with redshift of the scaling laws through the factors and A,. For 
instance, the Lboi - T relation (Eq.|9]l contains the cosmological term £ A^^. One definition of evolution (hereafter referred to as 
Evi) is any change with redshift in the scaling laws not accounted for by the and A, factors. A second definition of evolution 
(hereafter referred to as EV2) accounts for any change with respect to the local observed relations, including and A-. 

In the case of Evi the approach to follow is to correct the X-ray luminosity for the E^ and A, factors and then to introduce a factor 
(1 -I- z)'^- to account for any remaining effect with redshift. The Lboi - T relation then becomes: £'7'(A,/A,=o)"'^^Lboi oc T"'(l + z)'^'. 
With this notation in the self-similar scenario A, = and no Evi is present. 

For Ev2 generally a factor (1 -Hz)'* is used which includes also the E, and A, factors. For the Lboi - T relation one would write: 
Lboi °^ T"'{1 + z)^. In an Einstein-de Sitter Universe the cosmological factor E^ (EqlT]) and the overdensity A^ (EqjSl are equal to 
(1 + z)'^ and 1, respectively. The self-similar scaling relation would be: Lboi - CT"(1 + z)'^, and therefore A = 1.5. We would talk 
of evolution according to the definition of Ev2, but no evolution according to the definition of Evi . 

In a ACDM model with - 0.3, the dependence of E,aI^^ is no longer a power law of (1-i-z) but it can be approximated with that 
law in a limited range of redshifts by the redshift dependent exponent A : 



log 



(£-Io)(a,,;o) 



1/2 



log(l +z) 



where £,=0 = 1 (10) 



In the case of evolution Evi according to the redshift interval one can calculate the range of predictions for A on the basis of the 
self-similarity model and evaluate any possible evolution with respect to self-similarity. For example for a redshift range 0.3-hI.O 
the self-similar model prediction for A is in the range 1 -H 1.1. 



2.3. Breaking the Self-Similar Model 

One way to break the scaling laws that predict Lboi is to have non gravitational energy injected into the intracluster medium 
before or during the cluster formation. An abundant literature exists on this subject (for a review see e.g. rV^l2005al:lTbzzill2006l 
and reference therein). The candidates for the energy excess are currently feedback from star formation processes and feedback 
from nuclear activity in the cluster galaxies. Another physical process that could break the simple self-similar scaling is the removal 
of low-entropy gas via cooling (e.g. Bryan 2000; Voit & Bryan 2001). 

For example. IVoitI (12005 ah gives the following Lboi - T relations (his equations 81, 82, 83) for three different models: 
Lboi T^ HE,)-' 
Lboic^T2-'(£z)-'f(z)-' 
Lboi oc T\E,)-h(zr^ 

where t(z) is the cosmic time at redshift z. The above relations were obtained by using a cluster radius that encloses a mean density 
which is a fixed factor ( A - 200) with resp ect to the critical density. This alternative choice of a fixed overdensity is often found 
in the literature (see e.g. lEvrard et al. 1'2002'). In this case the factor A, is constant and the redshift dependence in the self-similar 
scaling relations (Eq.sO is only through the factor E^. 

The first model is a pure pre-heating case, in which the minimum entropy is assumed to be independent of both cluster mass and 
redshift. The other two are cooling models. These models give a slope of the Lboi - T relation significantly steeper than the self- 
similar model. They also predict different dependences on redshift (approximately oc (1 -h z) " (1 + z)" "*, and almost constant in 
the redshift range 0-1). The relations differ from the self-similar ones in their dependence on T and on z. Therefore, in principle, 
any observed change in the scaling laws allows to constrain different scenarios of heating schemes and cooling efficiency. 
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3. The cluster sample 

In this work we use Chandra archival data of eighteen distant (0.25 < z < 1.01) clusters. Details of the X-ray ana lysis are given in 
a companion paper dBranchesi et al.ll2007al Paper I). This same sample was also used bv lBranchesi et al.] (l2007bl) to check for the 
presence of overdensity of point sources in the inner region of clusters of galaxies. 

Since th e Lboi - T relationsh ip at both low and high redshift exhibits a large intrinsic scatter, mostly due to the strong cooling flow 
clusters (iFabian et al.ll 19941) . the presence of cooling cores has to be considered. To reduce the intrinsic scatter, which precludes 
the accurate determination of the exact shape and of any possible evolut i on of the Lboi - T relation, one can e xclude from the 
sample the clusters with strong cooling cores (e.g. Arnaud & Evrar d 1999; Ettor i et al.ll20 04; Mau ghan et"an ^2Q06). Alternatively, 
since cooling flow clusters are similar to non-cooling flow clusters except for the inner small cooling flow region, an approach 
to follow is to excise t he cluster central cores and replace the m with an extrapolation of the isothermal y6- model before deriving 
cluster parameters (e.g. lMarkevitchll998HVikhlinin et al ]|2002 l). There are four clusters in our sample which ' Vikhhnin et akl (|2002|) 
indicate as "possible coohng core" clusters. We have kept these clusters in the sample without applying any correction to T or 
to Lboi since according to their radial profile derived by us these cooling cores, if present at all, are rath er mild. A fifth cluster , 
ZWCL 1454.8+2233 presents a very thermally complex ICM (see notes on individual clusters in Paper I. iBranchesi et al.ll2007al) 
and was thus excluded from the analysis. From now on the sample of 17 clusters is referred to as '17 cluster sample'. 
In order to increase the statistics we supplemented ou r data w ith the sample of Vikhlinin et al.l ( l2002h (from now on VI02), 
lEttori et akl (12004.) (from now on ET04) and of Maughan et al.l ( 12006) (form now on MA06) after appropriate corrections. For 
clusters in common we use our own measurements. For clusters not in our sample but in common with the other three samples, we 
preferred the data of MA06 when available, since they use the most recent calibrations of Chandra. We also used some of the MA06 
clusters which were observed with XMM-Newton. For clusters common to VI02 and ET04 sample s, we use the VI02 clu sters since 
their data are in better agreement with ours after corrections for the different radii are made (see IBranchesi et ani2007a . Paper I). 
This choice was made "a posteriori" since the analysis evaluation of the combined data, given in Sect. |4] shows a worse value 
when using ET04 data, even if the final results do not change independently of the sample used. The final working sample is thus 
composed by 17 clusters from our sample, 8 clusters from VI02, 10 clusters from MA06 and 4 clusters from ET04 for a total of 39 
clusters. From now on the sample of 39 clusters is referred to as 'combined sample'. 

Since we want to compare our observational results with theoretical models, the choice of the radius within which integrated cluster 
properties are measured becomes very important, notably when a large range of redshits is involved. Several authors have chosen 
a constant radius which has the great advantage of simplicity. A constant radius would allow the analysis of the same fractional 
volume of clusters if clusters have all the same physical radius. However, in the self-similar models the virial radius is a function 
of redshift and temperature, thus a constant physical radius does not include a constant fraction of the cluster, and introduces a bias 
with redshift in the computed cluster luminosities. Since th e self-similar models pr edict the outer boundary of the virialized part of 
clusters in terms of a density contrast Av(z) (see for example fervan & Normanll998 ^. we consider appropriate to calculate the sample 
luminosities using the radius where the mean enclosed density is a factor A- above the critical density of the Universe (Ettori et al] 
1200 4; Maughan et al. 2006). In this way a fixed fraction of virial radius is used, independent of the cluster redshift, with no redshift 
and temperature biases in the luminosities. This choice allows us to use in a straightforward way the dependence of cluster properties 
on redshift as reviewed earlier (see Sect. 12. 11 12. 3t . In the present analysis we use for the overdensity A~ - 500 x Av(z)/(187r^) as in 
ET04. Note that the A-^o value is ^ 280, while the overdensity in a flat A = cosmology would be 500. For each cluster the value 
of R500 is calculated according to Eq.|6] All luminosities are extrapolated to Rsoo using an isothermal y6-profile. 

The uncertainties on R500 and L500 are calculated using the propagation of the measurement uncertainties on temperature and 
luminosities (which are the dominant contribution) and of the uncertainties related to the /3-model. Since core radius, Rc, and /3 are 
correlated quantities (surface-brightness model with large values of Rc and p are similar to those with small values), the uncertainties 
on them are not independent. Thus the errors on Rc were ignored in the calculation of the uncertainties (see MA06). 
It has to be noted that there may be some inconsistencies when comparing the high-z combined sample with the local samples due 
to a different choice of the radius used for the luminosity estimates. Markevitch ( 1998) adopted a radius of 1.4 Mpc, independent 
of temperature and /3, which is however close to R^m for average temperatures and p values due to the low redshift of their clusters. 
An evaluation of the ef fects due to different choices of the radius resulted in negligible differences as compared to the errors. 
lArnaud & EvrardI d 19991) do not specify any radius. We guess that they used a value encompassing all the visible X-ray emission. 
Even in this case we expect that systematic effects are negligible. 

The temperature and the Lsoo luminosities for all the clusters used to analyze the Lboi - T relation are listed in Table [T] It is worth 
to note that no dependence on redshift is present for either the luminosity or temperature of the sample clusters. The columns in 
Table [T]contain the following information: 

- Column 1: Cluster name. Asterisks indicate "possible cooling core" clusters as classified by VI02. However, no correction was 
applied to the four objects in the '17 cluster sample' since we checked that the cooling core, if any, is negligible. 

- Column 2: Spectroscopic redshift 

- Column 3: Cluster temperature in keV. For VI02, ET04 and MA06 the temperature is taken from their papers. 

- Column 4-5: Core radius in kpc and For VI02, ET04 and MA06 these parameters are taken from their papers. 

- Column 6: Radius in kpc which corresponds to a redshift-dependent contrast A, = 500 x Av(z)/(187r^). For each cluster the 
radius R500 is computed according to Eq. |6] 

- Column 7: Luminosity within R500 expressed in 10'*"'erg s"'. For the '17 cluster sample' the estimated luminosities (see 
feranchesi et al.ll2b07ai Paper I) are extrapolated to R500 using a yS-model. In the same way the VI02 luminosities (estimated 
within the fixed radius of 1.4 h^y Mpc) and the MA06 luminosities (estimated within the radius R200 corresponding to an 
overdensity A, = 200 x Av(z)/Av(0)) were extrapolated to R500 
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- Column 8: Original sample 



Table 1. Combined sample: cluster parameters 



Cluster name 


z 


T 


Rc 




R500 


L500 


Sample 


Abell2125 


0.246 


Q 4 + 0.2 

- «.? 


182+12 


54+0 02 


888+28 


2 13+0 08 


'17 rlnstpr' 


MS 1008.1-1224 


0.302 


6 + 
4 4 + 0.5 

75 ^^^^ 

- 0.7 


165+9 


64+0 02 


1241+42 


11 89+0 37 




ZW CL 0024.0+1652* 


0.394 


128+10 


67+0 03 


1003+58 


4 38+0 25 




MS 1621 5+2640 


0.426 


227±17 


65+0 03 


1255+84 


12 98+0 69 

1. • ^ Kj -i-\J -yj Zf 




RXJ 1701.3+6414* 


0.453 


4.5 + "-^ 


15±2 


0.41+0.01 


765+33 


7.14+0.40 




CL 1641+4001 


0.464 


5.1 


151+18 


0.77+0.06 


1092+91 


3.02+0.25 




V 1524.6+0957 


0.516 


9.4 + 


302+27 


0.80+0.05 


1029+72 


7.10+0.52 




MS 045 1.6-0305 


0.539 


270+8 


0.90+0.02 


1505+51 


52.89+2.75 




V 1121+2327 


0.562 


4 5 ; 8J 

^•■^ - 0.4 


437+58 


1.19+0.18 


1108+116 


5.73+0.43 




MS 2053.7-0449* 


0.583 


4 3 + !!-^ 


115+12 


0.64+0.03 


839+49 


5.81+0.55 




V 1221+4918 


0.700 


' " - 0.7 


272+20 


0.76+0.04 


1037+64 


13.22+0.79 




MS 1137.5+6625 


0.782 


6.2 + 
3.7 + 11 


116+6 


0.71+0.02 


910+39 


13.95+0.53 




RDCSJ 1317+2911' 


0.805 


61 + 16 


0.52+0.04 


592+87 


1.35+0.40 




RDCSJ 1350+6007 


0.805 


4.1 + 


261+43 


0.70+0.07 


688+77 


5.33+0.78 




RXJ 1716.4+6708 


0.813 


6.5 + °| 

- 0.8 


119+11 


0.66+0.03 


882+61 


13.35+1.08 




MS 1054.4-0321 


0.830 


8.3 + 


520+32 


1.38+0.11 


1339+89 


35.73+1.79 




WARPJ 1415.1+3612 


1.013 


6.2 + I* 

— U. / 


68+7 


0.60+0.02 


722+43 


12.27+0.88 




CL 1416+4446* 


0.400 


3.7 + ^li 


152 


0.66+0.04 


906+47 


4.14+0.29 


VI02 


MS 0302+1658* 


0.424 


3.6 1^ 

- 0.5 


154 


0.74+0.09 


929+88 


5.15+0.37 




RASS 1347-114* 


0.451 


14.1 + 

- 0^ 

5.3 + "-^ 
9.9 + 11 


84 


0.65+0.01 


1706+56 


135.01+9.45 




3C295* 


0.460 


103 


0.67+0.03 


1051+55 


8.04+0.56 




MS 0016+1609 


0.541 


274 


0.74+0.02 


1396+41 


57.38+4.03 




CL 0848+4456 


0.574 


2.7 + ti 

• ' - 0.3 


78 


0.58+0.04 


636+42 


1.29+0.09 




RDCS 0910+5422 


1.100 


3 5 + " 

^- ' - 1.0 


111 


0.72+0.18 


550+92 


2.66+0.21 




RDCS 0849+4452 


1.261 


119 


0.85+0.33 


628+144 


2.95+0.22 




RXJ 0542-4100 


0.634 


7 Q + 1.0 

' - 0.8 


132+17 


0.51+0.03 


982+63 


12.15+1.36 


ET04 


RXJ 2302+0844 


0.734 


6-6 ' oi 


96+12 


0.55+0.03 


865+65 


5.45+0.17 




RXJ 1252-2927 


1.235 


5 2 + 0" 

•'■^ - 07 


77+13 


0.53+0.03 


532+40 


5.99+1.10 




RXJ 0848+4453 


1.270 


9 Q + 0.8 
^•^ - 0.8 


163+70 


0.97+0.43 


499+115 


1.04+0.73 




CIJ 0046.3+8530 


0.62 


4 4 + U.5 

^•^ - 04 


137+28 


0.60+0.06 


792+56 


3.82+0.22 


MA06 


CU 1342.9+2828 


0.71 


2 7 + 0.5 


172+26 


0.70+0.06 


724+55 


3.34+0.18 




CIJ 1113.1-2615 


0.73 


4 7 + 8:t 

^- ' - 0.7 


106+12 


0.67+0.04 


801+74 


3.73+0.32 




CIJ 1103.6+3555 


0.78 


6 + 


141+18 


0.58+0.11 


807+93 


4.62+0.30 




C1J0152.7-1357N 


0.83 


249+48 


0.73+0.10 


820+93 


10.68+1.42 




C1J0152.7-1357S 


0.83 


- 1.0 


123+24 


0.66+0.07 


744+93 


6.42+0.96 




CIJ 1559.1+6353 


0.85 


4.1 


67 +30 


0.59+0.09 


646+107 


2.47+0.31 




CU 1008.7+5342 


0.87 


3.6! 4 

10.6! 1:1 


170+43 


0.68+0.09 


622+79 


3.57+0.35 




CIJ 1226.9+3332 


0.89 


113+8 


0.66+0.02 


1070+58 


42.51+1.00 




CU 1429.0+4241 
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4. The Lboi-T relationship for high-z clusters: the analysis method 

In order to examine how cluster luminosity scales with temperature and whether the Lboi-T relation depends on redshift, we followed 
four different steps in our analysis. 

We first assumed the following evolution law (Evi): 

Lbol,44 = CT«(l+z)^ (11) 

where Lboi,44 is the bolometric luminosity in units of 10'*'* erg s * and Tg - T(keV)/6. The parameters C, a and A are determined 
by fit which takes into account the uncertainties on the luminosity and temperature (ctl, (Tt)- The;^^^ is written as: 



2 _ V [logLboi,44,i - logC - glogTe^i - A log(l + a)]^ 

^ ~ / 1 2,92 ^ ' 



i=l 



where N is the number of clusters, ep is equal to cr-xjiyi:, In 10) and eL is equal to (TL/(Lln 10). The denominator represents the 
propagation of the uncertainties on the measured quantities. The minimum is found by varying log C, A and a in a grid of values. 




Fig. 1. {Left panel) The Lboi - T relation for the sample of 17 distant clusters. The solid line indicates the best fit to the data with A 
set to 0. Dotted and dashed lines refer to the nearby cluster samples. (Right panel) Same as above but for the 'combined sample' of 
39 clusters which includes data from MA06, VI02 and ET04 as explained in the text. 



a) The fit is performed setting A equal to zero. In this way we analyze the Lboi - T relation independently of the cluster redshift. 

b) The three parameters A, a and log C are used as free parameters in the fit. 

c) In order to constrain the evolution in the Lboi - T relation with respect to nearby galaxy clusters, an approach similar to VI02 
and ET04 is adopted. The fit i s performed by fixi ng i n Ea.[T2]a and log C to t he best fit values (a, log C) obtained from two 
local samples, namely those bv lMarkevitchI d 1 998h and lArnaud & Evrardl ( ll999h . lMarkevitchl ( IT99 8) analyzes 30 nearby clusters 
with measured ASCA temperature s and j^O&4numinositie s. His cluster luminosities are computed within a constant radius (1 .4 
Mpc for //() = 70 Km s ' Mpc"^). lArnaud & EvrardI (Il999t) analyze 24 nearby clusters with X-ray measurements from GINGA, 
A SCA and Einstein. B oth samples have a median re dshift of 0.05. In our cosmology the Lboi - T slope and normahzation found 
bv lMarkevitchI (Il998l) and lArnaud & EvrardI d 19991) convert respectively to: 

a = 2.64 + 0.27 and log C = 0.80 + 0.04 

a = 2.88 + 0.15 and log C = 0.77 +0.03 (13) 



As a last step the luminosities are scaled by the cosmological factor E. '(A,/A,=o) '''^ as predicted by the self-similar model (see 
e.g. ET04 and MA06). The data are fitted using the following evolution law (Evi): 

£-'(A,/4=o)"'^'Lboi,44 = CT^d + zf^ (14) 
d) The minimum is found either for A, equal to zero and leaving A, , a and log C as free parameters. 

With this approach the redshift dependence expected by the self-similar scenario is removed 'a priori' and thus normalization and 
slope values consistent with those from the local relations should be obtained if the self-similar predictions on the evolution are 
correct. 



5. The Lboi-T relationship for high-z clusters: results 

5. 1 . Comparison witii local samples 

The left panel of Fig. [T] shows the Lboi - T relationship for the '17 cluster sample' while the right panel shows the same plot for 
the 'combined s ample' of 39 clusters. Ov erplotted are the best-fit relationships found for the nearby clusters by Markevitch (1998|) 
(dotted line) and lArnaud & EvrardI d 19991) (dashed line). A displacement of the distant Lboi - T relation from the local ones is clearly 
evident, in agreement with the results of VI02, ET04 and MA06. High redshift clusters are on average systematically more luminous 
than the local ones by a factor ^ 2. 

To quantify the differences between the high and low-z clusters the Lboi - T fitting procedure, described in steps a) and b), has been 
applied to the '17 cluster sample' and to the 'combined sample'. The results of all the fits are listed in Table|2] Setting A=0 (first 
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Fig. 2. Ratio of the observed luminosity to the 'expected luminosity' from the local relations (see Sect. 15. Il l versus (1 + z) for 
the 'combined sam ple'. {Left panel) The expected luminosity (L^Jj'^j^^ly^jj.jj) is the luminosity obtained from th e local relation of 
iMarkevitchI 119981) using the observed cluster temperature. (Right panel) Same as above but for the local relation of Arna ud & EvrardI 
(1999) . In both panels the solid line indicates the fitted evolution law according to step c) of the Sect.|4] The evolution laws are 
clearly not appropriate. The error bar at z=0 indicates the uncertainty on the normalization of the local sample best-fit law. 



and third line in Table |2l), the best-fit values for the slope a are somewhat higher but still consistent within the errors with the local 
ones (see Eq.s[T3]l and significantly steeper than the self-similar predictions (a - 2). The best-fit values for log C diSei from the 
local ones at more than 5 cr. 

When the Ev2 evolution factor (1 + z)'^ is introduced in the fit (second and fourth line in Tabled, the;^'^ slightly improves but a 
and C do not change significantly suggesting that the dependence on redshift must be weak. In fact the Ev2 parameter A does not 
differ significantly from zero. Therefore there is no evidence for a change in the normalization of the Lboi - T relation as a function 
of redshift in the range covered by the 'combined sample' (0.25 < z < 1.3) or, in other terms, no Ev2 is present. 

T he data were t hen fit ted ( step c)) by fixing the slop e and normalization to the best-fit relations obtained for the low-z cluster samples 
of lMarkevitchI ( 119981) and lArnaud & EvrardI (1 19991) . The results are summarized in Table [3] A significant positive evolution (of Ev2 
type) with respect to the local relations is required: A varies between 1.2 and 1.35 for the 'combined sample'. A similar positive 
evolution is found by other authors (e.g. VI02, Kotov & Vikhlinin 2005, Lumb et al. 2004, and MA06i3- However, the best-fits are 
significantly worsened with respect to the previous step b) analysis. In particular the;^'^ of the fit for the 'combined sample' has a 
probability lower than 0.1% to be acceptable. This evolution, which would be consistent with the self-similar model predictions 
(A X 0.9 1.2 for the redshift range 0.25 H- 1.3 with respect to redshift 0), is not found in the high-z sample, but only when the 
entire redshift range (0 < z < 1.3) is considered, that is using the local samples as well. 

In order to emphasize this result, we present in Fig. |2]the ratios between the obser ved cluste r luminosity and the 'expected luminos- 
ity' derived using the observed high-z cluster temperature in the local relations bv iMarkevitch (. 1 998.) and Arnaud & EvrardC 19991) 
as a function of (l-Hz). For each cluster of the 'combined sample' the 'expected luminosity' is calculated as: 

T expected T^T^ M 

^Markevitch.Ainaud ~ ^ 6 ^^-'^ 

where a and C are fixed to the best-fit values obtained bv' Markevitc h' ( 1 998*) and lArnaud & E vrard'('l999') (see Eq.slTTb. respectively, 
and T is the measure of the high-z cluster temperature. Clearly the evolution law fitted according to step c) (see Table |3]l and 
indicated by the solid line is not a good representation of the data. 



' VI02 and Kotov & Vikhliniiil ( 120051) found A = 1 .5 and A =1.8, respectively, somewhat larger than our value. This is largely acco unted for by 
their different choice for the radius (1.4 Mpc) within which Lboi is computed. Using their same radius and comparing our sample with iMarkevitch! 
( !l998i) , we find A =1.6 in good agreement with them. 
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T (keV) T (keV) 

Fig. 3. The £7' (A./A-=o) ''^^ Lboi - T relation for the '17 cluster sample' (left panel) and for the combined 39 cluster sample (right 
panel). The luminosity of each cluster is scaled by the cosmological factor £7' (A-/A,=o) The solid line indicates the best-fit to 
the data obtained fixing A, - 0. The dotted and dashed lines indicate the local relations. 



5.2. Comparison with fhe self-similar model 

In Table m we list the results of the fits obtained by scaling the luminosities by the cosmological factor £7' (Ar/A2=o) as 
described in step d), and fixing A, - 0. The best fit values of parameters a and C are found to be very similar to those of the low 
redshift clusters. This result is illustrated in Fig. |3] where the high redshift cluster data points, corrected for the cosmological factor, 
are randomly distributed on the low-z relations (see also ET04 and MA06). However the;^'^ value in Table |4] for the 'combined 
sample' shows that the fit is not good (probability lower than 1% to be acceptable). This is emphasized in Figure|4]that shows the 
ratios between the luminosity scaled by the cosmological factor £7' (A,/A,=o) '^^ and the luminosity corresponding to the cluster 
temperature according to the local relations as a function of (1+z) for the 'combined sample'. If the evolution were consistent with 
the predictions of the self-similar models, then the data points should be randomly distributed around unity, independent of the 
redshift. Figure |4] shows instead an excess of clusters above unity forz < 0.6 and a large excess below unity for z > 0.6. 

When Ay is let free to vary we find a strong negative Evi evolution with the normalization factor C increasing by a factor 2 (see 
TableHli which brings the data points to shift back, away from the local relations. This is the same effect mentioned earlier. In order 
to bring the high-z clusters data points in agreement with the local relation a redshift dependent coiTection must be applied to the 
data. This effect is not present in the redshift range explored by the combined sample of 39 clusters. 

We are led to think that the dependence on redshift of the Lboi - T relation in the full range < z < 1.3 cannot be described by 
either the factor E~ (A. / A-^ q)'^^ or by a p ower law factor ( 1 + z)'* . These contradictions between data and evolution models (Ev i and 
EV2) were already noted bvl VoitI (l2005bh who used a collection of published data. Voit pointed out that the luminosity evolution can 
not be monotonic with redshift. He suggested that the self-similar evolution may be plausible at low redshift, but cannot continue 
to arbitrarily high redshift. This might be consistent with radiative cooling and feedback from galaxy formation being increasingly 
more important at high redshift. 

Unless there are systematic observational effects between the low-z and high-z cluster data (for instance a systematic temperature 
difference of 25% or an appropriate combination of systematic differences in luminosity and temperature) our results require a 
strong evolution, similar or stronger than the self-similar, from z = 0toz<0.3 followed by a null evolution at higher redshift. 
Systematic diffe rences in both lu minosity and temperature between past and current generation of X-ray detectors are not reported 
in the literature. I Vikhhnin et al.l (i200Z) specifically verified that no bias is present between the temperatures measured by Chandra 
and by ASCA. 

An additional effect could be the presence of the Malmquist bias in flux-limited samples (Sect. 5.2 in lMaughanll2007al) . However, 
this effect is very difficult to quantify with the existing heterogeneous archival cluster samples. 
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Fig. 4. (Le/f panel) Plot of the ratio bet'ween the luminosity scaled by the cosmological factor Ez^ (A./Aj=o) '^^ and the 'expected 
luminosity' corresponding to the cluster temperature according to the local relation of [M arkevitch (19981) versus (1 + z) for the 
'combined sample'. (Right panel) Same as above but for the local relation of lArnaud & EvrardI (119991) . In both panels the error bar 
at z=0 indicates the uncertainty on the normalization of the local sample best-fit law. 



5.3. Comparison witii additional Cinandra data 

A very recent 'work bv ' Maughan et alJ ( 2007b|) analyzes a sample of 1 15 galaxy clusters at 0.1 < z < 1.3 observed 'with Chandra 



It is interesting for us to use the lMaughan et al] (l2007bh sample since they have a large number of nearby objects (z < 0.4) which 
are missing in our sample. From their 1 15 clusters we selected 56 clusters as follows. Twenty-five objects are in common with our 
'combined sample' clusters and, according to us, are not contaminated by the central "cooling core". These 25 clusters were used 
to verify the presence of systematic difference between ours and Maug han measurements. O ur estimates of the cluster temperature 
and luminosity are on average 15% and 10% higher with respect to Maughan et al.l (l2007bl) . We do not discuss here the possible 
reasons for these systematic differences, but just take them into account in the analysis. The remaining clusters were selected from 
the z < 0.4 sample and are those clusters for whic h the ratio between the luminosities with and without the central cluster region 
is less than 1 .7 (see columns 4 and 8 of Table 3 in iMaughan et al. I l2007bh . In such a way we exclude clusters with strong cooling 
cores. Figure |5] is the same as Fig. |2] (left panel) but for the Maughan's cluster measurements. A comparison between Fig. |5] and 
Fig. |2] shows that the two distributions are very similaj^ and confirms and extends to lower redshift our results. In fact the z > 0.25 
clusters have on average a flat distribution higher than LMaughan/LMa^i^evftch ~ ^ which is the local cluster value. At lower redshift (z 
< 0.25) the region pertaining to the local samples begins to being suddenly populated (8 clusters). The reason of this change in the 
distribution is not at all clear and not accounted for to our knowledge by any model. 



5.4. Comparison witin non-gravitational models 



Figure 



6] shows the comparison between the high-z 'combined sample' and the three non-gravitational models described 



(l2005ah where the effects of preheathing and radiative cooling are included in the simple self-similar scenario. To properly compare 
our data to the Voit models the luminosities were extrapolated to a radius corresponding to a fixed overdensity (A = 200) with 
respect to the critical density as required by the models. The three panels in Fig|6]give the ratio of the extrapolated luminosity to the 
T", where a is the slope of the Lboi - T in each model versus (1h-z). The dashed lines represent the fit of the models to the data, the 
dot-dashed lines represent the case when the models are normalized to the local samples at z=0. 

It is worth noting that if only the 'combined sample' (0.25 < z < 1.3) is considered the three models (dashed lines in Fig.|6|l are in 
agreement with the cluster data. The model that better describes the redshift evolution of the Lboi - T relation is the third model ( 
Lboi K T^(E^)^^t(z)^^) shown in the bottom panel of Fig.|6] The corresponding minimum has a probability of about 50% to be 
acceptable (x^.^/d.o.f. = 37.3/38) against a probability of about 5% for the other models (model 1 xlji„/d.o.f. - 51.8/38, model 2 



^ The systematic differences in luminosity and temperature mentioned above convert into a higher normalization (about 25%) in the y-axis of 
the Maughan distribution with respect to ours. 
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Fig. 5. Ratio of the luminosity measured bv lMaughan et alJ (l2007bh to the 'expected luminosity' derived using the lMaughan et alJ 
(I2()07bl) measured temperature in the local relations by iMarkevitchI (Il998l) for the 56 clusters (see text for details). Filled circles 
indicate the 25 objects in common with our 'combined sample'. Open circles are objects with z < 0.4 a nd with a ratio betwee n the 
luminosities with and without the central cluster region less than 1 .7 (see columns 4 and 8 of Table 3 in lMaughan et al ] l2007bh . 

xfjji^/d.o.f. = 57.3/38). The third model has also the same temperature slope (a = 3) as the fit of a power law (La=200 = CTg(l +z)^) 
to the data (solid line in the bottom panel of Fig|6]l. 

When we link the high-z 'combined sample' to the local samples (see dot-dashed lines in Fig.|6l) none of these model is acceptable. 
Namely none of the non-gravitational models seems to be able to account for the redshift dependence in the entire redshift range 
<z < 1.3. 



6. Summary and Conclusions 

In this paper we have presented a re-visitation of the observational relation between X-ray luminosity and temperature for high-z 
clusters of galaxies. To this purpose we have constructed the Lboi-T relation for a sample of 17 distant clusters ( 0.25 < z < 1.01) 
selecte d from the Chandra archive. A detailed analysis of the X-ray data are given in a companion paper (Paper I bv lBranchesi et al.) 
I2007ah . The statistics was increased by adding to our sample 22 clusters taken from the literature to form a final sample of 39 
objects in the redshift range 0.25 < z < 1.3, observed either by Chandra or XMM-Newton. The analysis of the larger sample allows 
a more stringent determination of the high-z Lboi-T relation for a comparison with the local one and with the theoretical models. 
We emphasize the importance to extrapolate the luminosities to a fixed fraction of the virial radius when observational results 
are compared to theoretical predictions. The present analysis indicates that the apparently different results on the Lboi-T found by 
several authors can be in some cases explained by the different methods used, for instance a different choice of the cluster radius, 
a different definition of evolution (Evi or EV2) and/or a different way to compare the observational results with the self-similar 
predictions. 

Our analysis of the Lboi-T relation of high redshift clusters has revealed a significant evolution with respect to the local Universe. 
The evolution goes in the direction of high-z clusters being brighter than the local ones by a factor w 2 at any given temperature. 
The slope of the relationship is found to be steeper th an expected from the self-similar model pr edictions and steeper, even though 
still compatible within the errors, than the local slope (lMarkevitchlll998tlArnaud & Evrardll 19991) . 

For the analysis of the Lboi-T evolution with redshift, different approaches were adopted (see Section|4]l: 

- the data corresponding to the distant clusters in the redshift range 0.25 < z < 1.3 were fitted considering the evolution factor 

(l+z)-*; 

- a second fit is performed with the evolution factor ( 1 +z)'^ and fixing the slope and normalization to the best-fit parameters obtained 
for the well studied low-z cluster relations; 

- the luminosities are corrected for the cosmological factor predicted by the self-similar model E^' (Az/Az=o) '''^ before the fit is 
performed. 

The following results are found: 

i) The high-z 'combined sample' shows no significant changes with z in the explored range 0.3 < z < 1.3. The fit of the data with a 
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Fig. 6. Ratio between the La=2()0 and T", where a is the temperature dependence of the Lboi - T relation found by IVoitI ( l2005ah . 
versus (1 + z). The dashed lines indicate the redshift dependence of the non-gravitational models, considering only the 'combined 
sample'. The dot-dashed lines represents the non-gravitational models when they are normalized to the local samples. The solid 
line in the bottom panel represents the fit obtained following step b) described in Sect. |4]and using La=200 (A = -0.55 + 0.66, 
a = 3.00 + 0.22, log C = 1.20 + 0.13, ;t'^jj.^/(i.o./. = 1.01). The error bar at z=0 in each panel indicates the normalization of the local 
samples. 



power law (1 + z)'^ yields for A a value not significantly different from zero. 
ii) The self-similar model is not consistent with the data. 

Hi) When the high-z 'combined sample' is linked to the local samples no acceptable fit is obtained either using a power law (1 +z)'^ 
or the self-similar scaling law. 

The data from a recent paper bv'Mau ghan et al.l ( l2007bl) confirm our results down to z ^ 0.25. In addition, we note in the range 0. 1 
< z < 0.25 the appearance of a number of clusters with the same properties as the local clusters. 

The non-gravitational models bv lwi(l2005al) are in good agreement with the data when the redshift range 0.3 < z < 1.3 is consid- 
ered. When the entire range < z < 1.3 is taken into account none of these models is acceptable. 

Summarizing we find that the dependence of the L^-T relation on redshift in the full redshift range explored (0 < z < 1.3) cannot 
be described by the factor Ej (Az/Az=o)'^^ of the self-similar model, or by any power law of the form (1 -n z)'*. It seems that a strong 
evolution is required by the L^-T relation from z = to z ~ 0.3 followed by a much weaker, if any, evolution at higher redshift. This 
weaker evolution is c ompatible with an increasing importance at high redshift of non-gravitational effects in the structure formation 
process ( IVoitll2005bl) . 

According to us the first priority for a correct analysis of the evolution of the Lboi-T relation is to assemble a more statistically 
significant sample of nearby (z < 0.25) clusters. A handful of clusters not contaminated by cooling cores and at z < 0.1 could be 
enough to shed light on this issue. This should be kept in mind for the future generation satellites which should have a larger field 
of view and the same or better resolution and sensitivity than the present X-ray telescopes. 

An additional effect to take into account since it could mimic evolution is the presence of the Malmquist bias. With the available 
archival samples it is very difficult to keep under control and factor out this bias. It would be extremely useful to observe with 
Chandra or XMM-Newton complete samples of nearby and distant clusters all coming from the same survey with well defined 
selection criteria in order to be able to control and eventually correct for the bias. 
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sample 


A 




a 


logC 


Xl.jd-o.f. 


17 cluster sample 





3.14 


+ 0.34 

- 0.29 


1 09 + " 

'^ ■"^ - 0.04 


12.18/15 




-0.08!":^^ 


3.15 


+ 0.57 
- 0.40 


1 10 + 

- 0.18 


12.15/14 


combined sample 





3.00 


+ 0.19 
- 0.17 


1 06 + " °3 

- 0.03 


39.29/37 




-0.48 


3.03 


+ 0.25 
-0 20 


1 15 +"15 


37.18/36 



Table 3. Best fit results of Lboi - T relation - Luminosities extrapolated to R500 (step c) 



sample 


A 


Xljd.o.f. 


Reference sample 


17 cluster sample 


1.40!!J;|^ 

1.58 


22.39/16 
23.26/16 


Markevitch ^998) 
Amaud&Evrard(1999) 


combined sample 


1 20 + " 

- 0.08 

1 + 0.08 

- 0.10 


73.57/38 
68.35/38 


Markevitch (1998) 
Arnaud(&Evrard(1999) 



The slope {a) and the normalization (log C) are fixed to the best-fit values obtained for the low-z cluster sample indicated in the last column. 
Table 4. Best fit results of Er' (AJA;=o)"'^^ Uoi - T relation (step d) 



sample 


-4; 


a 


logC 


xL/d.o.f. 


17 cluster sample 




- 1.28 ! 


9 QQ + 0.33 
- 0.28 

T 17 + 0.55 
3.11 _ 


Q go + 0.04 
- 0.04 


17.99/15 
12.18/14 


combined sample 




-1.80!™ 


3 00 + " " 
3 05 + " 23 


0-87 

1.20!!;-" 


61.09/37 
37.37/36 



The luminosities extrapolated to R500 correspond to a redshift-dependent contrast of 500, and were scaled by the cosmological factor 
£■7' (A,/Aj=o) ''^ as predicted by the self-similar model. 
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